Strongly correlated quantum dots in weak confinement potentials and magnetic fields 
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We explore a strongly correlated quantum dot in the presence of a weak confinement potential 
and a weak magnetic field. Our exact diagonalization studies show that the ground state property 
of such a quantum dot is rather sensitive to the magnetic field and the strength of the confinement 
potential. We have determined rich phase diagrams of these quantum dots. Some experimental 
consequences of the obtained phase diagrams are discussed. 
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Fabrication of quantum dots has reached such an ad- 
vanced state that the shape of the confinement poten- 
tial and the number of electrons in them can be tuned 
precisely 0. Ground states of isolated quantum dots 
fabricated from 2D electron systems can display vari- 
ous strong correlation effects. Quantum Hall related 
ground states have been explored in the strong magnetic 
field limit 0,13, and Hund's rule for shell structures for 
a two-dimensional harmonic potential has been investi- 
gated 0. Recently transport through a quantum dot 
has attracted much attention, especially Kondo related 
physics Wigner crystal states and strong electron 
correlation in quantum dots have been explored through 
various methods |a,l3 

In this paper we explore possible ground states of an 
isolated quantum dot made of dilute 2D electron sys- 
tems in the presence of a weak magnetic field. Re- 
cently unusual behavior suggestive of a metal-insulator 
transition has been reported in a variety of dilute two- 
dimensional electron and hole systems 0, where the di- 
mensionless parameter R = E e - e /E s ^> 1 (here S e -e 
is the electron-electron interaction energy and E s is the 
characteristic single particle energy). Strong electron 
correlation effects are expected to be play an impor- 
tant role in these systems. Quantum dots fabricated 
from dilute 2D electron systems are also expected to ex- 
hibit strong electron correlation effects. Here we focus on 
how the ground state spin depends on the applied mag- 
netic field and the strength of the confinement potential. 
To describe properly the delicate competition|9( between 
different spins states in strongly correlated regime non- 
perturbative methods^, El [Dl 0] are required. We 
adopt exact diagonalization methods 13. [Tol fllj since var- 
ious mean field theoretical methods [l3j are applicable 
only in the regime R < 1. 

We model a quantum dot as follows. Electrons move on 
a 2D plane under the influence of a parabolic confinement 
potential and a weak magnetic field applied perpendic- 
ular to the plane. In parabolic quantum dots, the ratio 
R can be characterized by {e 2 /Ka)/Ml where a is the 
typical lateral size of the dot, k the dielectric constant, 
and f2 the confining frequency of the harmonic potential. 
We consider up to six electrons in the parameter regime 
R ~ 5 — 16. We find that the ground state of quantum 



dots for these values of R is rather sensitive to the mag- 
netic field and the strength of the confinement potential. 
Our investigation suggests that this is a direct conse- 
quence of strong electron correlation: it originates from 
the existence of nearly degenerate quantum dot eigen- 
states. We obtain rich phase diagrams of these quantum 
dots and discuss experimental consequences. 

Electrons of the dot is confined by a harmonic poten- 
tial V(r) = ^m*U 2 r 2 , where m* is the effective electron 
mass. A magnetic field B is applied along the z-axis 
through a vector potential in a symmetric gauge. We in- 
clude also the Zeeman splitting, g/.iB[meV] = 0.026S[T]. 
We take the Hamiltonian of N electron dot to be H = 
52i=i Hi + Hi n t, where the single particle Hamiltonian is 
modeled by 

H i = ^ + l m *" 2 r? - \uj c \i ■ B - gnBS,, i% (1) 
and where the many-body interaction term is given by 
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Here k — 12.4 and ui c are the dielectric constant of 
GaAs semiconductor and cyclotron frequency. The 
eigenfunctionspTL ITU of the single particle Hamiltonian 
are 



4>nUz (r) 
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where a 2 = H/(2muj), uj 2 = ft 2 + \io 2 c , \s z are spin func- 



tions, and 



R n i(x) 
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Note that, as expected, the dependence on 9 indicates 
that electrons rotate clockwise when a magnetic field ap- 
plied along the z-direction. However, it is convenient to 
define the z-component of single particle angular momen- 
tum as \ z = — 14 so that the eigenstates with positive 
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angular momenta have lower energy. The eigenenergies 
are 

e n i Sz = hw(2n + \l\ + 1) - ^hu c £ - g\iS z B. (5) 

These eigenstate wavefunctions are labeled by quantum 
numbers of orbital states and the z-component of angular 
momentum (n, I): Note that for a given n = 0, 1, 2, ... the 
quantum number I can take all possible integer values. It 
is useful to study properties of zero magnetic held single 
particle states. The eigenenergies are 

e° nl = m(2n + \l\ + l). (6) 

These single particle energy levels e° ; can be degenerate, 
for example, (n, /) and (n, — I). These degeneracies are 
such that they lead to the magic numbers of harmonic 
potential 2, 6, 12, etc. In the presence of a weak magnetic 
held along the z-axis all degeneracies will be lifted. 

The many-body ground state \P of H is expanded in 
terms of Slater determinant states 

* = £)C a * a , (7) 

a 

where 

Here a = {niliai, unInCn} and labels a Slater de- 
terminant. Each creation operator c^ ls creates an elec- 
tron in an eigenstate of zero magnetic field single particle 
Hamiltonian. The energy of ty a is denoted by E ai and 
it is the sum of energies of single particle states (\> n is z 
Even in the presence of a magnetic field we use the same 
many-body basis states as in zero magnetic field, with 
replacing £1 by ui. We found numerically that it is conve- 
nient to use these zero field particle states in constructing 
Slater determinants. The many-body Hilbert space can 
be divided into regions with different quantum numbers 
(S Z ,L Z ): In each region we construct the many-body ba- 
sis states which have the well-dehned z-componcnt of to- 
tal angular momentum quantum number L z = Y^iLi h 
and the z-component of total spin S z . All possible eigen- 
states in the Hilbert space (S Z ,L Z ) are calculated, and 
the minimum energy state is found. 

For numerical diagonalizaton of the Hamiltonian ma- 
trix, only Slater determinant basis states with energies 
E a < qhuj are included, where q = 20— 25. Typically, we 
have included about 3,000-14,000 basis states. We have 
checked that for these values q the energy differences be- 
tween different states have converged. For example, when 
N = 6, B = and Ml — 0.3meV the energy differences 
between the ground state (0, 0) and one of the compet- 
ing state (2, 0) are 0.0371, 0.0332, 0.0328 meV for q=21, 
23, 25, respectively. Results for N = 3 and 4 are in 
good agreement with the results of Mikhailov 7]- For 
example, for Ml = 0.118565 (corresponding to A = 10 
in Ref. [7]), we have E {1/2A) /hn = 17.6279(iV = 3) and 
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FIG. 1: Phase diagram for N = 3. Each state is denoted 
by the quantum numbers (S Z ,L Z ). For strongly correlated 
regime (for small h£l), the ground states are very sensitive to 
the magnetic field. 
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FIG. 2: Same as in Fig.l, but N = 4. 
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FIG. 3: Same as in Fig.l, but N = 5. 
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FIG. 4: Same as in Fig.l, but N = 6. 

£(i i0 )/Ml = 31.4120(iV = 4), perfectly consistent with 
the numbers in the reference. Note that when Ml gets 
too small more basis states need to be included and the 
numerical work becomes unwieldy. These results are not 
included in the phase diagrams. 

Figures 1,2,3 and 4 show phase diagrams for N — 
3, 4, 5, and 6 in the weak potential and magnetic field 
regime. Each phase is labeled by the total spin and 
angular momentum quantum numbers (S Z ,L Z ). These 
phase diagrams are constructed schematically from the 
calculated energies of different (S Z ,L Z ) states at various 
values of {Ml,B). We took B = 0, 0.01, 0.19, 0.2 T 
for all values of N shown in the figures, and took Ml = 
0.1, 0.2, ... 1.0 mcV for N = 3,4,5 while Ml = 0.2, 0.3, 
0.5, 0.7, 1.0 meV for N — 6. We remark on some of 
the qualitative aspects of these phase diagrams. (i)Note 
that some the phase regions are divided into subregions 
separated by horizontal phase boundaries. Also we no- 
tice that some of the phase boundaries are linear func- 
tions of B. (ii)These figures show that, as the strength of 
the potential or magnetic field changes the ground state 
quantum number L z changes sensitively. As the mag- 
netic field B increases the value of L z tend to increase. 
On the other hand as the strength of the confinement po- 
tential Ml increases the ground state L z tend to decrease. 
(iii)Ground state energies obtained for L z = 0,3,4, and 
6 are, respectively, 37.28, 37.58, 37.88, and 38.05 meV 
(These values are for the parameters N — 6, B — 0, and 
Ml = 1.0 mcV). These results indicate that the ground 
states are nearly degenerate as a result of strong electron 
correlation. 

We now discuss why some of the phase boundaries in 
Figures 1, 2, 3, and 4 arc linear functions of B. It is 
useful to separate the total energy into 

B±ot — ^in ~t~ Ediam "h Bp ara -\- E S pi n -\- E{ n t , (9) 

where Ekin =< e n i- > i s the kinetic energy, 
Ediam =< J2i l m * UJ c r i > is tne diamagnetic energy, 
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FIG. 5: Here we compare the differences E pa ra{S z ,2, L z , 2) — 
E P ara(S Z: i, i z ,i) and E (5 z ,i, L z ,i) - E (S z ,2, L Z)2 ) for states 
sharing the same phase boundary. Two solid lines repre- 
sent the magnetic field dependence of AE para (AL z = 4) = 

E P ara(0,0)-E P ara(l,4:) and Ae' 4() = e' (1, 4) - E' (0, 0) , and 
the other two dotted lines represent AE para {AL z = 2) = 
E para (l,4) - E para (0,6) and AE' 64 = E (0,6) - E'(l,4). 
When these pair of lines cross the ground state changes. Here 
iV = 6, hQ, = 0.7meV, and E spin is ignored. Note that the 
differences in E (S Z ,L Z ) are almost constant as a function of 
magnetic field while the difference in E par a{S z , L z ) are signif- 
icant. 



Epara = —\u c L z is the paramagnetic energy, E int is 
the interaction energy, and E spin = —g[iBS z is the 
Zeeman energy. Note that Ekin oc ft. We now com- 
pare Ep ara with the rest of the ground state energy, 
E ' (S Z ,L Z ) = E kin + E diam + E int , ignoring E spini which 
is very small compared to the others. For many states 
found in the phase diagrams the values of E (S Z ,L Z ) 
arc nearly the same. This property reflects the deli- 
cate competition between different states, which is an 
indication that the quantum dot is in the strongly cor- 
related regime. However, since Ep ara depends signifi- 
cantly on B, we expect ground states to change sensi- 
tively as a function of B. This is explained in detail 
in Figure 5. Phase boundaries in the parameter space 
(HSl, B) is roughly determined by the sum of kinetic and 
paramagnetic energies: Consider the phase boundary be- 
tween ground states [S Zt i,L z ,i) and (5 Z ,2, L z ^), then, 
since the values of E {S Z ^,L Z ^) are nearly the same, 
we expect bVl\ — ^uj c L Zt i « b£l2 — \i^ c L z ^, which gives 
fii — £1 2 ~ (L z ,i — L z _2) (here b is the proportionality 
factor in the kinetic energy). This result indicates that 
some phase boundaries are approximately a linear func- 
tion of B. Also the slope of such a phase boundary is 
proportional to the difference in the angular momentum 
L Zi \ — L Zt 2- These results are consistent with many phase 
boundaries in Figures 1, 2, 3, and 4, although there are 
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FIG. 6: Radial density profiles of several different ground- 
states with S z = 0. Here Ml = 0.7meV and N = 6. For 
weak magnetic fields, radius of dot is roughly 5a, where 
a — ^h/2mu>. We notice that the profiles deviate slightly 
from each other. From these profiles we can estimate the 
density of 2D electron gas which these dots are made of to be 
about 9 x 10 9 cm" 2 . 



few exceptions. 

We now discuss some issues of experimental relevance. 
First we estimate the density of 2D electron gas which 
these dots can be made of. Figure shows the ra- 
dial electron density profile of dots for N = 6 when 
Ml = 0.7meV. The radius of the dot is roughly 5a for 
weak fields considered in this paper. We have verified by 
integrating the density profile that the total number of 
electrons is N. For m* — 0.067m, where m is the bare 
electron mass, the average electron density is approxi- 
mately 9 x 10 9 cm -2 . We suggest that our obtained phase 
diagrams may be explored experimentally by measuring 
the positions of cusps in the magnetic field dependence of 
chemical potential hn = En — -Ejv-i: Since ground state 
level crossings in the N — 1 and N electron systems lead 
respectively to positive and negative cusps, phase bound- 
aries will show up as cusps in the chemical potential. We 
plot the magnetic field dependence of the chemical po- 
tentials /i4, /its, and He, in Figured We observe several 
cusps in Fig. [71 where Ml = 0.7 meV is used. From our 
phase diagrams we can conclude that for smaller values 
of Ml many more cusps will show up in the magnetic 
field dependence of the chemical potential. Conductance 
peak spacings, given as the difference between the gate 



issue concerns that in real dots there may be some de- 
viations from the perfect circular symmetry of the har- 
monic potential, and consequently degenerate single par- 
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FIG. 7: Plot of the chemical potential /ijv = En — En -i as a 
function of the magnetic field at hfl = 0.7 meV for N = 4, 5, 
and 6. The positions of magnetic fields, where a ground state 
level crossing in the N — 1 and N electron systems takes place, 
are indicated by A and V- The quantum states of the dots 
are denoted by (S Z ,L Z ) (The other cusp- like structures are 
artifacts of discrete data points). 



ticle states of the harmonic potential at B = may be 
absent. However, even in a harmonic potential electron- 
electron interactions will lift this degeneracy. Further- 
more, the main conclusion of our work that the ground 
states change sensitively with B or Ml is not expected to 
be different since strong electron correlation is responsi- 
ble for the sensitivity of ground states to the magnetic 
field and the strength of the confinement potential. 

In this paper we explored quantum dots in a strongly 
correlated regime R 3> 1. Our exact diagonalization re- 
sults show that the ground states of such quantum dots 
are rather sensitive to the magnetic field and strength of 
the confinement potential. We have predicted rich phase 
diagrams. Experimentally these phase diagrams may be 
explored by measuring the sensitive magnetic field depen- 
dence of the energy to add one electron to a dot, or by 
measuring the positions of conductance peak oscillations. 
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